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Abstract. The existence of cosmic rays of energies exceeding 10 20 eV is one of the mysteries of 
high energy astrophysics. The spectrum and the high energy to which it extends rule out almost 
all suggested source models. The challenges posed by observations to models for the origin of 
high energy cosmic rays are reviewed, and the implications of recent new experimental results are 
discussed. Large area high energy cosmic ray detectors and large volume high energy neutrino 
detectors currently under construction may resolve the high energy cosmic ray puzzle, and shed light 
on the identity and physics of the most powerful accelerators in the universe. 
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1. Introduction 

Fig. 1 presents a schematic description of the spectrum and composition of cosmic rays ob- 
served at Earth [1]. At low energies, ~ 1 GeV per particle, the flux is dominated by protons. 
The average particle mass increases with energy and the composition becomes dominated 
by heavy nuclei at ~ 10 6 GeV per particle, where the spectrum also becomes steeper. At 
still higher energy, ~ 10 10 GeV, the spectrum and composition change again: the spectrum 
becomes harder (flatter) and there is strong evidence that the composition becomes lighter, 
likley dominated by protons. It therefore appears that a new source of cosmic rays comes 
to dominate above 10 19 eV. Since heavy nuclei of energies < 10 18 eV are confined by the 
Galactic magnetic field, it is believed that cosmic rays of energy < 10 19 eV are of Galactic 
origin. This view is supported by the (small but statistically significant) enhancement, ob- 
served below 10 19 eV, of cosmic ray flux from the direction of the Galactic plane. Protons 
of energy > 10 19 eV are not confined by the Galactic magnetic field, and the isotropic 
distribution of cosmic rays above 10 19 eV therefore suggests that the flux is dominated at 
these energies by an extra-Galactic source of protons. 

The origin of the highest energy, > 10 19 eV, cosmic rays (UHECRs) is a mystery. As 
explained in §2., the high energies observed rule out almost all candidate sources. The 
situation is further complicated by the interaction of high energy protons with microwave 
background photons. As explained in §3., this interaction limits the propagation of pro- 
tons of energy > 10 20 eV to < 100 Mpc, and there are no exceptionally bright sources 
that may be suspected as UHECR sources within such a distance from Earth. More over, 
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Figure 1. A schematic description of the differential spectrum and of the composition 
of cosmic rays observed at Earth. 

it is not clear whether or not the expected "GZK suppression" [2] of UHECR flux above 
~ 5 x 10 19 eV, due to interaction with microwave background photons, is observed. These 
difficulties have led to the suggestion that modifications of the basic laws of physics are 
required in order to account for the existence of UHECRs. Such suggestions are discussed 
in detail in a separate contribution to these proceedings [3], The present contribution fo- 
cuses on the discussion of possible solutions to the UHECR puzzle, which do not invoke 
modifications to the basic laws of physics. In §4. we discuss the gamma-ray burst (GRB) 
model for UHECR production and some of its predictions, which may be tested with giga- 
ton neutrino telescopes. A more detailed discussion of the model and its predictions for 
planned UHECR and neutrino detectors may be found in [4]. Some general comments on 
the role that neutrino telescopes may play in resolving the UHECR puzzle are given in §5.. 
Our main conclusions are summarized in §6.. 



2. Acceleration: General considerations and candidate sources 

A detailed discussion of particle acceleration is beyond scope of this talk. However, the 
essence of the challenge can be understood using the following simple argument. Consider 
an astrophysical source driving a flow of magnetized plasma, with characteristic magnetic 
field strength B and velocity v. Imagine now a conducting wire encircling the source at 
radius R, as illustrated in fig. 2. The potential generated by the moving plasma is given by 
the time derivative of the magnetic flux $ and is therefore given by V — (3BR where j3 = 
v jc. A proton which is allowed to be accelerated by this potential drop would reach energy 
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Figure 2. Potential drop generated by an outflow of magnetized plasma. 

outflow, where T = (1 — /3 2 ) -1 / 2 ^> 1. In this case, the proton is allowed to be accelerated 
only over a fraction of the radius R, comparable to R/T. To see this, one must realize that 
as the plasma expands, its magnetic field decreases, so the time available for acceleration 
corresponds, say, to the time of expansion from R to 2R. In the observer frame this time is 
R/c, while in the plasma rest frame it is R/Tc. Thus, a proton moving with the magnetized 
plasma can be accelerated over a transverse distance ~ R/T. This sets a lower limit to the 
product of the magnetic field and source size, which is required to allow acceleration to 
Ep, 

BR > TE p /e(3. (1) 

Figure 3 shows B and R for various objects. GRBs are the only type of sources which 
safely satisfies the basic requirement of eq. 1 . Active galactic nuclei (AGN) appear in two 
places: as "AGNs", where acceleration is considered near the central source (massive black 
hole) which drives the outflow, and as "Radio Galaxies", where acceleration is considered 
in regions of the flow far away from the central object, where interaction with surrounding 
gas leads to radio emission. Acceleration to high energy near the central source is impos- 
sible due to the high density of photons, which leads to rapid energy loss of protons by 
photo-production of pions. Acceleration in galaxy clusters is impossible since the flows 
there are characterized by /3 < 10~ 25 . 

Eq. 1 also sets a lower limit to the rate L at which energy should be generated by the 
source. The magnetic field carries with it an energy density B 2 /8tt, and the flow therefore 
carries with it an energy flux > vB 2 /8tt (some energy is carried also as plasma kinetic 
energy), which implies L > vR 2 B 2 . Using eq. 1 we find 

Only two types of sources are known to satisfy this requirement. The brightest steady 
sources are active galactic nuclei (AGN). For them V is typically between 3 and 10, im- 
plying L > 10 47 crg/s. This is marginally satisfied by the brightest AGN. The brightest 
transient sources are GRBs. For these sources T ~ 10 25 implying L > 10 50 5 erg/s, 
which is satisfied since the typical observed MeV-photon luminosity of these sources is 
i 7 - 10 52 erg/s. 
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Figure 3. "Hillas plot" [5]. The magnetic field strength and size of different objects. 
In order to allow proton acceleration to E p , sources must satisfy eq. 1. 

It was recognized early on ( [5] and references therein) that while highly magnetized 
neutron stars may lie above the required line in the "Hillas plot" (fig. 3), it is hard to uti- 
lize the potential drop in their electro-magnetic winds for proton acceleration to ultra-high 
energy. The mechanisms recently proposed for acceleration in "magnetars" ( [6] and refer- 
ences therein) also face serious difficulties. In [6], for example, the electro-magnetic wind 
must penetrate through a supernova envelope shell without losing energy to acceleration of 
the shell and without "contaminating" the wind with baryons. The mechanism by which 
such penetration may be achieved is unclear. 

3. Propagation: The GZK suppression 

3.1 Interaction with microwave background photons 

As illustrated in Fig. 4, high energy protons may interact with cosmic microwave back- 
ground photons to produce pions. In each interaction of this type, the proton loses a frac- 
tion ~ m^/irip of its energy. The threshold energy requirement, E V E 1 > m p m T c 4 where 
E p and E 7 are the proton and photon energies respectively, implies that protons of energy 
E p > 10 20 eV may interact with almost all of the T — 2.7° K background photons, while 
protons of lower energy may interact only with the tail of the Planck distribution. Thus, 
the energy loss distance, Xe(E p ), drops rapidly with energy in the range of 0.5 x 10 20 eV 
to 3 x 10 20 eV (see Fig. 4). 
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Figure 4. Schematic description of the strong reduction of the energy loss distance, 
\e, of high energy protons, due to pion production interactions with microwave back- 
ground photons. E20 is the proton energy in units of 10 20 eV. 

The rapid drop of Xe(E p ) with increasing E p has two major consequences. First, it ef- 
fectively limits the distance out to which sources may contribute to the UHECR flux above 
10 20 eV to less than ~ 100 Mpc, a small distance on cosmological scale. In particular, 
there is no AGN bright enough to satisfy the constraint of Eq. 2 with this distance. This 
appears to rule AGN as the sources of UHECRs. This difficulty may be avoided by assum- 
ing the existence of "dead AGN" [7], massive black holes in nearby galaxies that produce 
UHECRs without producing almost any light. The difficulty here is that even if a small 
fraction of the energy output is converted to electron, rather than to proton, acceleration, 
the resulting radiation field may inhibit proton acceleration (through energy loss to pion 
production). The implications of the distance limit to GRBs are discussed in detail in [4]. 

The second implication is that for any homogeneous distribution of sources of UHECRs, 
the drop of Xe{E p ) with increasing E p will strongly suppress the flux above ~ 0.5 x 
10 20 eV: At lower energies, sources out to distances ~ 1 Gpc, comparable to the size of 
the observable Universe, may contribute to the flux, while at higher energies only local, 
nearby sources contribute to the flux. This suppression of the flux, usually named the 
"GZK suppression," is generic to all models where UHECRs are protons produced by a 
homogeneous distribution of extra-Galactic sources. The question of whether or not such 
suppression is present in the data is therefore extremely important. The analysis below, 
aimed at answering this question, follows that of ref. [8], where interested readers may 
find more details. 



3.2 Analysis of current data 

In order to address the question of the presence or absence of a GZK suppression, we 
must adopt a model for the intrinsic spectrum produced by the UHECR sources (which is 
subsequently modified by the GZK effect). We assume that extra-galactic protons in the 
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energy range of 10 19 eV to 10 21 eV are produced by cosmologically-distributed sources at 
a rate and spectrum given by 

« °- 6 x lo44er S M P C " 3 yr~VW. (3) 

Ct-C/p 

An energy spectrum similar to the assumed dN/dE p oc E~ 2 has been observed for both 
non-relativistic [1] and relativistic [9] shocks. It is believed to to be due to Fermi accelera- 
tion in collisionless shocks [1,9], although a first principles understanding of the process is 
not yet available. The normalization of 0.6 x 10 44 erg Mpc -3 yr _1 is chosen to account for 
the observed flux. (j){z) accounts for redshift evolution (4>(z = 0) = 1). Motivated by the 
GRB model (§ 4.), we have assumed that <fi(z) follows the evolution of the star-formation 
rate. The spectrum above 10 19 eV is only weakly dependent on <p(z) since proton en- 
ergy loss limits their propagation distance. For the heavy nuclei component dominating at 
lower, < 10 19 eV, energy we take the Fly's Eye experimental fit [10], 

f oc E-^. (4) 
dE 




E [eV] 



Figure 5. Model versus data. The solid curve shows the energy spectrum derived 
from the two-component model discussed in § 3.2. The dashed curve shows the ex- 
tra-Galactic component contribution. The "No GZK" curve is an extrapolation of the 
E~ 2 - 75 energy spectrum derived for the energy range of 6 x 10 18 eV to 4 x 10 19 eV 
[1]). Data taken from [11]. AGASA's energy scale was chosen. 

Figure 5 demonstrates that model predictions are in good agreement with the data of all 
experiments in the energy range 10 19 eV to 10 20 eV 1 . Above 10 20 eV, the Fly's Eye, HiRes 



As explained in detail in [8], the various experiments are consistent with each other when system- 
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and Yakutsk experiments are in agreement with each other and the model. However, the 
eight AGASA events with energies greater than 10 20 eV disagree with the prediction of 
the cosmological model. The Fly's Eye, Yakutsk and HiRes experiments have a combined 
exposure three times that of the AGASA experiment. The exposures above 10 20 eV are, 
in units of 10 3 km 2 - yr - sr: AGASA (1.3), Fly's Eye (0.9), Yakutsk (0.9), and HiRes 
(2.2). Together, Fly's Eye, Yakutsk, and Hi-Res observe a total of 6 events above 10 20 eV 
(4 events if the Fly's Eye energy scale is chosen). 

How significant is the flux suppression observed by Fly's Eye, HiRes and Yakutsk? The 
differential energy spectrum observed by the various experiments at the energy range of 4 x 
10 17 eV to 4 x 10 19 eV can be fitted by a broken power-law, where the shallower component 
dominating above ~ 6 x 10 18 eV satisfies J oc E~ 2 - 75±0 - 2 . The number of events observed 
beyond 10 20 eV by Fly's Eye, HiRes and Yakutsk is at a > 5<r deficit relative to the number 
expected from extrapolation to high energy of the low-energy distribution. Adopting the 
steepest allowed slope, J oc E~ 2 - 95 , the observed number of events is at a > 3.7a deficit. 

What is the reason for the discrepancy between AGASA and the other experiments above 
10 20 eV? Our analysis shows that the discrepancy is statistically significant. Thus, the 
discrepancy is most likely due to systematic errors in the estimates of event energies. The 
Auger detector currently under construction [ 1 2] is likely to resolve this issue by combining 
two different detection methods (ground array and nitrogen fluorescence). 

4. GRBs, UHECRs and neutrinos 

Gamma-ray bursts (GRBs) are short, typically tens of seconds long, flashes of gamma-rays, 
carrying most of their energy in > 1 MeV photons. The detection in the past few years 
of "afterglows", delayed X-ray, optical and radio emission from GRB sources, proved that 
the sources lie at cosmological distances, and provided strong support for the scenario of 
GRB production described in fig. 6 [13]. The energy source is believed to be rapid mass 
accretion on a newly formed solar-mass black hole (or, possibly, neutron star). Recent 
observations suggest that the formation of the central compact object is associated with 
type Ib/c supernovae [14]. 

The energy release drives an ultra-relativistic, T ~ 10 2 5 , plasma outflow. At a large 
distance from the central black-hole, internal collisionless shocks within the wind, which 
arise due to variability in the wind emitted from the central "engine", accelerate electrons 
to high energy. Synchrotron emission from these shock accelerated electrons is believed to 
account for the observed 7-rays. At still larger distances, the wind impacts on surrounding 
medium. Here too, the collisionless shock driven into the ambient gas accelerates electrons, 
leading to synchrotron emission which accounts for the "afterglow". 



atic errors in the absolute energy scale of the events are taken into account. The relative systematic 
shifts in absolute energy calibration between Fly' Eye and the other experiments, required to bring 
into agreement the fluxes measured at 10 19 eV by the different experiments, are {-11%, +7.5%, 
-19%} for {AGASA, HiRes, Yakutsk}. All shifts are well within the published systematic errors. 
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Figure 6. The fireball scenario of GRB production. 



Wind property 

u B /u e 

r 

dN/dE 



MeV 7-rays 

> 0.1 

> 300 

dN e /dE e oc E~ 2 



10 20 eV protons 

> 0.02 

> 100 

dN p /dE p oc E- 2 



Table 1. Constraints on GRB wind parameters from photon observations and from the 
requirement for proton acceleration to 10 20 eV. 



4.1 The association ofGRBs and UHECRs 

GRBs were suggested to be UHECR sources in [15,16]. The GRB-UHECR association 
was based in [15] on two major arguments: (i) The constraints imposed on the relativistic 
wind by requiring that it would produce the observed 7-rays are remarkably similar to those 
imposed by the requirement that the wind would allow proton acceleration to 10 20 eV; (ii) 
The rate (per unit volume) at which energy is generated by GRBs in 7-rays is similar to the 
rate at which energy should be generated in high energy protons in order to account for the 
observed UHECR flux. 

The constraints on wind parameters are summarized in table 1 (for a pedagogical dis- 
cussion see [4]). 7-ray observations require the ratio of magnetic field and electron energy 
densities, ub/u c , to exceed 0.1 in order to produce electron synchrotron emission in the 
MeV range, V > 300 to avoid pair-production by high energy photon which would make 
the wind opaque, and dN e /dE e oc E~ 2 to account for the synchrotron spectrum. Accel- 
eration of protons to 10 20 eV require ub/u c > 0.02 to satisfy Eq. 1 (We have seen that 
this equation implies a lower limit to the luminosity carried by magnetic field, eq. 2, and 
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since the observed photon luminosity represents the wind luminosity carried by electrons, 
a lower limit is inferred for the ratio of magnetic and electron luminosity); T > 100 to 
avoid synchrotron losses of protons on a time scale shorter than the acceleration time, and 
dNp/dEp oc Ep 2 to reproduce the observed UHECR spectrum (see eq. 3). 

The fact that similar constraints are obtained from independent arguments, suggests that 
GRBs and UHECRs are associated. The association is further supported by the similarity 
of energy generation rates. We have shown in § 3. that the energy generation rate of protons 
inferred from UHECR observations is given by eq. 3. The energy generation rate of GRBs 
may be estimated as follows. Assuming that the GRB rate follows the star-formation rate, 
which is reasonable given the association with type Ib/c supernovae, the local (z = 0) 
GRB rate is ss 0.5/Gpc 3 yr [17]. The average MeV 7-ray energy release is ss 3 x 10 53 erg 
for GRBs with known redshift [18]. The rate of bursts for which redshift is obtained is 
smaller than the total observed rate by a factor of « 0.7 x 0.5 = 0.35. The factor 0.7 is due 
to the fact that the detection threshold of the BeppoSAX detector that allowed afterglow 
detection is higher (by a factor w 2) than the threshold of BATSE [19], the detector that 
provided the data based on which the total rate estimate was derived. The factor of 0.5 is 
due to the fact that the fraction of bursts for which optical afterglow was detected is <~ 0.5. 
Thus, the local (z = 0) energy generation rate of GRBs in MeV photons is 

£ 7 [MoV] > 0.35 x 3 x 10 53 erg x 0.5/Gpc 3 yr = 0.5 x 10 44 erg/Mpc 3 yr. 

(5) 

This rate is similar to the energy generation rate in UHECRs, eq. 3. 



4.2 Predictions: High energy neutrinos 

Protons accelerated in the fireball to high energy lose energy through photo-meson in- 
teraction with fireball photons. The decay of charged pions produced in this interaction 
results in the production of high energy neutrinos. The key relation is between the ob- 
served photon energy, E 7 , and the accelerated proton's energy, E p , at the threshold of the 
A-resonance. In the observer frame, 

£ 7 £ p = 0.2GcV 2 r 2 . (6) 

For T « 300 and E 1 = 1 MeV, we see that characteristic proton energies ~ 10 16 eV are 
required to produce pions. Since neutrinos produced by pion decay typically carry 5% of 
the proton energy, production of ~ 10 14 eV neutrinos is expected [20]. 

The fraction of energy lost by protons to pions, / w , is f n w 0.2 [4]. Assuming that GRBs 
generate the observed UHECRs, the expected GRB muon-neutrino flux may be estimated 
using eq. 3 [4], 

E l*» ~ -J- f -T E l(. dl *p/dE p )t H « 0.3 x lO-^GeVcm-V^r" 1 . (7) 
47T 4 U.z 

Here tn is the Hubble time and the factor 1/4 multiplying f n is due to the fact that in photo- 
production of pions charged and neutral pions are created with roughly equal probability, 
and when a charged pion decays roughly half its energy is carried by muon neutrinos. This 
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neutrino spectrum extends to ~ 10 16 eV, and suppressed at higher energy due to energy loss 
of pions and muons. Eq. 7 implies a detection rate of ~ 20 neutrino-induced muon events 
per year (over 4ir sr) in a cubic-km detector. Since GRB neutrino events are correlated 
both in time and in direction with gamma-rays, their detection is practically background 
free. 

High energy neutrinos may be produced also in other stages of fireball evolution. For a 
detailed discussion see [4,21] and references therein. 



5. Giga-ton neutrino detectors and UHECRs 

If UHECRs are protons of extra-Galactic origin then, regardless of the nature of their 
sources, their existence implies the existence of a flux of extra-Galactic high energy neutri- 




E [GeV] 



Figure 7. The upper bound imposed by UHECR observations on the extra-Galactic 
high energy muon neutrino intensity (lower-curve: no evolution of the energy pro- 
duction rate of eq. 3, upper curve: assuming evolution following star formation rate), 
compared with the atmospheric neutrino background and with the experimental upper 
bound of BAIKAL [23] and AMANDA [24]. The curve labelled "GZK" shows the 
intensity due to interaction with micro-wave background photons. Dashed curves show 
the expected sensitivity of O.lGton (AMANDA, ANTARES [25], NESTOR [26]) and 
lGton (IceCube [27], NEMO [28]) Cerenkov detectors, and of the balloon radio exper- 
iment ANITA [29]. 

nos. Such neutrinos are expected to be produced by the decay of charged pions produced 
in interactions of UHECRs with photons (or nucleons). UHECR observations set an upper 
limit to the intensity of high energy neutrinos produced in sources where the optical depth 
to nucleon interaction with photons is not high, i.e. where the average number of photo- 
production interactions the nucleon undergoes before escaping is not high [22]. This upper 
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limit is obtained by taking f v — 1 in eq. 7, i.e. by assuming that all the UHECR energy is 
converted to pions. The resulting bound, shown in fig. 7, applies to GRBs and to the ob- 
served jets of AGN, which satisfy the requirement of small photo-production optical depth. 
Sources with large optical depth may exist, and may therefore produce a flux exceeding the 
bound. However, we have no direct evidence for the existence of such sources. 

Fig. 7 demonstrates that the detection (by muon optical Cerenkov detectors) of the ex- 
pected fluxes in the energy range of 1 TeV to 1000 TeV requires cubic -km (i.e. Gton) scale 
detectors. The GZK flux (which "touches" the bound at ~ 10 19 eV because ultra-high 
energy protons produced at large distances lose all their energy to pion production in in- 
teraction with microwave background photons) requires still larger detectors (due to the 
lower number flux at higher energy). For a more general, yet concise, discussion of high 
energy neutrino telescopes and their science goals see chapters 2 and 3 of [30]. For more 
detailed review see [31]. 



6. Conclusions 

Detectors of high energy cosmic -ray and neutrinos currently under construction may allow 
to identify the sources of UHECRs. Such identification will provide only a partial reso- 
lution of the puzzle. A major challenge will remain in understanding the physics of the 
sources. GRBs and AGN are the most powerful astronomical objects, and are likely can- 
didates for the production of ultra-high energy protons and neutrinos. In both, the energy 
source is likely to be mass accretion onto a black hole, leading to relativistic outflows. The 
models describing these objects are largely phenomenological, and major open questions 
remain regarding the underlying physics. Data from the new experiments may allow to 
resolve some of these open questions. 

High energy neutrinos are expected to be produced in astrophysical sources by the decay 
of charged pions, which lead to the production of muon and electron neutrinos. However, 
if the atmospheric neutrino anomaly has the explanation it is usually given, oscillation to 
z/ T 's with mass ~ 0.1 eV [32], then one should detect equal numbers of v^s and ^ r 's. Up- 
going t's, rather than /it's, would be a distinctive signature of such oscillations. Since ^ T 's 
are not expected to be produced, looking for t's would be an "appearance experiment." 

Detection of neutrinos from GRBs could be used to test the simultaneity of neutrino and 
photon arrival to an accuracy of ~ Is, checking the assumption of special relativity that 
photons and neutrinos have the same limiting speed. These observations would also test 
the weak equivalence principle, according to which photons and neutrinos should suffer the 
same time delay as they pass through a gravitational potential. With 1 s accuracy, a burst 
at 1 Gpc would reveal a fractional difference in limiting speed of 10~ 17 , and a fractional 
difference in gravitational time delay of order 10~ 6 (considering the Galactic potential 
alone). Previous applications of these ideas to supernova 1987A (see [33] for review), 
yielded much weaker upper limits: of order 10~ 8 and 10~ 2 respectively. 
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